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Abstract: Image denoising and segmentation play a crucial role in computer graphics 
and computer vision. A good image-denoising method must effectively remove noise 
while preserving important boundaries. Various image-denoising techniques have 
been employed to remove noise, but complete elimination is often impossible. In this 
paper, we utilize Partial Differential Equation (PDE) and generalised cross-validation 
(GCV) within Adaptive Haar Wavelet Transform algorithms to effectively denoise an 
image, with the digital image serving as the input. After denoising, the image is 
segmented using the Histon-related fuzzy c-means algorithm (H-FCM), with the 
processed image serving as the output. The proposed method is tested on images 
exposed to varying levels of noise. The performance of image denoising and 
segmentation techniques is evaluated using metrics such as Peak Signal-to-Noise 
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Ratio (PSNR) of 77.42, Mean Squared Error (MSE) of 0.0011, and Structural 
Similarity Index (SSIM) of 0.7848. Additionally, segmentation performance is 
measured with a sensitivity of 99%, specificity of 98%, and an accuracy of 98%. The 
results demonstrate that the proposed methods outperform conventional approaches in 
these metrics. The implementation of the proposed methods is carried out on the 


MATLAB platform. 


DOI: 10.52756/ ijerr.2024.v44spl.005 


Introduction 

The transmission of multimedia substances such as 
digital images has recently been widely developed on the 
internet (Islam et al., 2016). In hyperspectral images 
(HSD), the eminence of the image is significantly affected 
by the occurrence of noise. So, denoising is carried out as 
a spatial preprocessing procedure which is used to 
develop the image eminence for data analysis (Aswathy 
et al., 2015). Image denoising is a common problem in 
image processing. As a result of several types of noise, 
digital images 


are often blurry. In digital image 


processing, noise is a quantity that can change the 
experimental value from its original value (Bhandari et 
al., 2016). The digital image is somewhat distorted by 
noise during the transmission or acquisition of the image. 
The foremost objective of image denoising is to evaluate 
the original image, such as the edges, textures, and fine 
details of its noise-degraded analysis. In many image 
processing systems, image denoising is a crucial 
preprocessing step, and the exploration of function is an 
essential challenge (Shaoping et al., 2016). Additionally, 
it is utilized to enhance the prominence of images to help 
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visual assessment. MRI denoising is a difficult problem 
in medical image processing because of the improbability 
of resources and the allocation of noise (Ai, Danni et al., 
2015; Gautam et al., 2022; Jain et al., 2023; Goel et al., 
2023; Himabindu et al., 2024). 

Gaussian is typically used for noise recognition. 
Magnetic resonance (MR) images with noise have been 
found to be in Rician (Phophalia et al., 2014, 2016). 
Additionally, several research works for image denoising 
have recently been predicted. However, the challenge is 
still seen in the latest studies. Usually, a linear 
representation like a Gaussian filter is used to diminish 
the noise (Malini et al., 2015). The flat area of the images 
is where these procedures perform best. One of the main 
advantages of the linear procedure is its relatively high 
computing efficiency. However, it is difficult for these 
procedures to reconstruct the entire original image and 
achieve the optimal denoising properties using only the 
local information of the observed image (Feng et al., 
2015). Additionally, non-linear representation is 
employed to maintain limits that are superior to those of 
the linear model (Jalab et al., 2015; Naresh and 
Peddakrishna, 2023). 

Here, several additional algorithms are also employed 
to eliminate noise. One approach that derives from 
thresholding is the multiresolution denoising algorithm. It 
is used to keep the residual as unaltered (hard- 
thresholding) or decrease it to the threshold value by 
setting all wavelet coefficients to zero in the chosen 
threshold value (Huang et al., 2015, 2016, 2014). Here, 
the Gaussian smoothing filter has become the standard. 
This type of local averaging procedure is to eliminate 
noise with structural information such as anatomical 
limitations. The actual noise has a huge impact on how 
the denoising process is presented. The impact of noise 
on several procedures, including image segmentation. 
Image segmentation is essential and significant for 
computer visualization, pattern identification, and 
medical image processing. The foremost intention of 
image segmentation is to divide an image mechanically 
or semi-automatically as non-overlapping regions. Every 
area is standardized by some uniqueness like intensity or 
texture (Feng and Bhandari, 2016). The anisotropic 
diffusion function is a characteristic of the partial 
differential equation used for noise reduction and the 
preservation of image features, including edges. The 
work introduces a tangent sigmoid diffusion coefficient 
and an adaptive threshold value, which results in 
accelerated convergence (Kollem, 2024). An approach to 
image denoising utilizing a diffusivity function-based 
partial differential equation that takes advantage of the 
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statistical features of the observed noisy images. This 
model employs a Quaternion Wavelet Transform, which 
generates the various coefficients of a noisy image. The 
soft threshold function is employed to enhance the 
generalized cross-validation function, which identifies the 
optimal threshold value. The ideal threshold value 
subsequently regulates the diffusion process through a 
novel diffusivity function. They present the fourth-order 
partial differential equation diffusivity function, a 
distinctive diffusion coefficient that surpasses previous 
methods in noise reduction and edge preservation 
(Kollem, 2023). 

This work aims to create an efficient approach for 
image denoising and segmentation that addresses the 
shortcomings of current strategies with respect to noise 
elimination and boundary preservation. The proposed 
approach utilizes Partial Differential Equations (PDE) 
and generalised cross-validation (GCV) in the adaptive 
Haar Wavelet Transform to markedly diminish noise 
while preserving the image's structural integrity. The 
Histon-related fuzzy c-means algorithm (H-FCM) is 
employed for segmentation to improve the precision of 
image segmentation, especially in noisy conditions. The 
primary objective is to achieve enhanced performance in 
denoising quality and segmentation accuracy, evidenced 
by superior performance measures including PSNR, 
SSIM, MSE, sensitivity, specificity, and accuracy relative 
to traditional approaches. 

Literature Survey 

This section discusses the image denoising and 
segmentation techniques that various authors have 
proposed. The approach suggested by Kollem et al. 
(2022) uses a Quaternion Wavelet Transform to create 
the various noisy image coefficients, an improved 
generalized cross-validation function to produce the ideal 
threshold value, and a new diffusivity function to regulate 
the diffusion process. This approach is not superior in 
terms of eliminating noise and maintaining sharp edges. 

An approach to denoising that is based on anisotropic 
diffusion processes has been proposed by Ally et al. 
(2021). To effectively preserve edges, contours, and 
textures, this approach’s divergence and regularization 
terms include a diffusion kernel that depends on the 
evolving image and its gradient magnitude. In MRI 
images, this technique causes distracting artifacts. 

Lotfi et al. (2022) have presented a method that 
contains a partial differential equation, numerical 
approaches, and the Perona-Malik equation to eliminate 
noise in an image. This method reduces the size of the 
original problem by decomposing it into a set of smaller 
problems, each of which is addressed differently. Being 
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unidirectional also eliminates the necessity for the 
intricate process of domain decomposition while solving 
the equations. The Perona-Malik methods produce the 
desired denoised results using meshless techniques. 

An approach to maintaining the edges of the image 
has been suggested by Rafsanjani et al. (2022). This 
approach combines efficient diffusion-based technology 
with residual feedback. The suggested method employs 
residual feedback and the localization of the noisy pixels, 
and it is carried out in two stages. There is a loss of detail 
in the image and problems with large gradients when 
using this technique. 

Abdellahoum et al. (2021) have come up with a 
method that uses a cooperative approach and improved 
fuzzy c-means. This cooperation is controlled by a multi- 
agent system that makes it possible to choose the best 
metaheuristics parameters automatically. The 
disadvantage, which results in a reduction in the 
segmentation accuracy, is the improper selection of 
clusters and their initialization. 

Kollem et al. (2020) developed a method that 
combines centroid optimization using the elephant 
herding approach with Atanassov's intuitionistic fuzzy set 
Histon-based fuzzy clustering method. This optimization 
strategy is used to choose clusters effectively. The 
instability in clustering is a result of the method's use of 
clan updating and separation operators. 

The fuzzy c-means clustering was proposed by Mishro 
et al. (2020) to segment brain MR images. This method's 
algorithm for segmenting brain tissue is a type-2 adaptive 
spatial fuzzy c means technique. The fixed value of the 
fuzzifier employed in fuzzy c-means, which introduces 
uncertainty in regulating the fuzziness of the extracted 
clusters, is a downside of this method. 

Zhang et al. (2019) demonstrated image segmentation. 
using a hybrid approach Using an evolutionary algorithm 
to find improved clustering solutions, this technique first 
generates a random population of initial clustering 
solutions at random then, at each iteration, the FCM 
technique is applied to each of these solutions to generate 
its segmentation result. 

Liu et al. (2016) established an image denoising 
process by means of developing the sparse depiction of 
boundaries and the multidirectional reduction. The 
augmentation of sparsity was accomplished by means of 
implementing directionlet transforms, which are created 
by the directions of boundaries. Here, the created 
directionlet transforms are carried out in dissimilar 
directions for every pixel. They acquired numerous 
The 


dissimilar estimations, which were the finest. 
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concluding denoised output was acquired by means of a 
weighted averaging of entity evaluation. 

Saeid Saryazdi et al. (2016) specify that the Perona- 
Malik (PM) dispersion coefficient was investigated by 
flux. According to the investigation, the PM flux of the 
section, where the gradient magnitude was superior to the 
smoothing threshold, is directly related to unwanted 
blurring consequences and border dislocation. The image 
was separated into three segments to tackle those 
problems that are derived from the gradient magnitude: a 
section where the gradient was lesser than the smoothing 
threshold, a section where the gradient was among the 
smoothing threshold and modulation position of flux, and 
a section where the gradient magnitude was superior to 
modulation position. They describe the circumstances 
that are measured in those three segments. After that, the 
conditions 
coefficient. 


are satisfied by a calculated dispersion 


Khan et al. (2016) presented an innovative impulse 
noise recognition algorithm that is derived from noise 
proportion assessment and a mixture of K-means 
clustering and Non-Local Means related filter (NEK- 
NLM). Luo-statistic has engaged a nonlocal means- 
related estimator. The uniqueness of the process is 
handling the establishment of the preprocessing step of 
noise proportion assessment before noise recognition, 
which is used to pick the appropriate limitation for the 
noise recognition algorithm (Xiaoyan et al., 2016). The 
nonlocal-means estimator was implemented in the noise 
filtering phase to re-establish the noisy pixels for their 
definite values. In actual world datasets, the design 
illustrates that the impulse noise is detached efficiently. 

Zhang et al. (2015) established the anticipated 
algorithm, which employs partial differential equations 
(PDE) and generalized cross-validation (GCV) theory in 
the tetrolet transform field. It is used to reduce noise in an 
image. The noisy image was analyzed using the Tetrolet 
transform, and the finest denoising threshold in the 
Tetrolet transform domain was determined using the 
GCV theory. Subsequently, the inverse tetrolet transform 
was employed to the customized tetrolet transform 
coefficients to acquire a preliminary denoising image. 
The PDE representation was employed to diminish the 
block consequence that occurs by means of tetrolet 
transform to decrease noise and maintain the border 
information (Wei Xie et al., 2016). The anticipated image 
denoising algorithm was evaluated by way of five 
comparable image denoising algorithms such as wavelet 
transform combined with PDE, contourlet transforms 
combined with PDE and curvelet transform combined 
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PDE, shearlet transform combined with PDE and tetrolet 
transform combined with PDE correspondingly. 
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Figure 1. Architecture of the Proposed PDE Threshold Function-Based Image Denoising 
and Histon FCM-Based Segmentation. 


Problem Definition 

A review of the research shown above has some of the 
problems with recent image denoising and segmentation 
methods as follows: 

# The newly developed diffusivity function in 
(Kollem et al., 2022) suppresses the image details. 

# In this approach, the divergence and regularization 
terms include a diffusion kernel that gives blurry images 
(Ally et al., 2024). 

# This method employs effective diffusion-based 
technology with residual feedback, although there is a 
loss of image detail and issues with large gradients when 
employing this method (Rafsanjani et al., 2022). 

# This method makes use of a cooperative approach 
and an improved fuzzy c-means algorithm, both of which 
result in an incorrect selection of clusters and 
initialization of those clusters (Abdellahoum et al., 2021). 

# This methodology combines  Atanassov's 
intuitionistic fuzzy set Histon-based fuzzy clustering 
algorithm with the elephant herding strategy. Border 
displacement occurs as a result of clustering instability 
(Kollem et al., 2020). 

# This technique uses a type-2 adaptive spatial fuzzy c 
means algorithm to segment brain tissue. One drawback 
of this approach is the fixed value of the fuzzifier used in 
fuzzy which introduces 
controlling the fuzziness of the retrieved clusters (Mishro 
et al., 2020). 


c-means, uncertainty in 
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Contributions 

The following are the contributions to address the 
limitations identified in literature survey. 

# The tetrolet transform is used as an application to 
extract the image 
coefficients. 


into low and_ high-frequency 

# Improved generalized cross-validation is used to 
obtain the optimal threshold value. 

# A new diffusion coefficient is used in the partial 
differential equation based on Improved generalized 
cross-validation for image denoising. 

# An improved Histon-based fuzzy c means clustering 
is proposed for image segmentation. 


Proposed Methodology 

In this part, we have addressed the implementation of 
the proposed method. The proposed approach for image 
denoising has the primary objective of removing noise 
while preserving boundaries' integrity. The primary 
objective of the proposed image segmentation is to 
prevent over-segmentation and noise sensitivity. The 
following is a description of Figure 1: The preliminary 
MRI images of the brain tumor, which were obtained 
from the reference resource known as BraTS2021(Baid et 
al., 2021). The tetrolet transform is used to decompose a 
noisy image into low- and high-frequency coefficients. A 
threshold approach is used to eliminate noisy image 
coefficients. In the tetrolet transform, the GCV method is 
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utilized to locate the finest denoising threshold value at 
each scale level. The proposed partial differential 
equation method produces a denoised image using an 
optimal threshold value. With the help of a Histon-related 
fuzzy c-means method, the denoised image is segmented 
(Binge Cui et al., 2016). As a result, we get segmented 
images as output. The sensitivity, specificity, and 
segmentation accuracy of PSNR, SSIM, MSE have all 
been used to calculate the performance of the proposed 
approaches. The proposed method is briefly discussed 
below: 

Tetrolet transform 

Using partial differential equations (PDEs) and 
generalized cross-validation (GCV) in the tetrolet 
transform, a novel image-denoising algorithm is expected 
to be presented in this article. The proposed algorithm's 
complexity is reduced by using the tetrolet transform's 
admirably sparse representation. Based on tetrominoes, 
the Tetrolet Transform is a novel adaptive wavelet 
transform in the Haar category. In general, the 
tetrominoes encompass four equal-sized squares, and 
each boundary is linked to another. 

Consider the following scenario: for each N*N image, 
if the value of N is even, then the image may be enclosed 
using only the five fundamental tetrominoes. At the 
beginning of the process, the image is divided into 4*4 
blocks at each stage of the tetrolet transform field. After 
that, using a technique called tetrominoes segmentation, 
which is determined by the uniqueness of the block's 
geometry, four sections are extracted. The input image is 


U =(uli, jJ)xj9 and where N=2’(JeN) . The 


equations for the low-frequency and high-frequency 
coefficients can be found below. 


sla 
u, =(u, li. i])?,-6 and 


wlijl= DY oxi.) h EF] 


i,j evi, 


() 


The low-frequency coefficients are calculated with 
Equation (1). The matrix of the tetrolet transform is 


denoted by the letter K (i, i). while the coefficient of 
the tetrolet transform [28] is denoted by the letter 77 . 


ae) 


a (u, [i. ee 


ufiil= DO ol bK(7) ks] 


(Jeu, 


and 
(2) 


The high-frequency coefficients are calculated with 
Equation (2). Where /=1,2,3 denotes the various 
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components that make up the high-frequency 
coefficients. 

The steps involved in the tetrolet transformation 
(Kollem et al., 2019) are outlined below. 
USC.) where 


Input: and 


N=2' (JEN) 

Step-1: During this process stage, the image is 
divided into 4*4 blocks. 

Step-2: Equations (1) through (2) are used in the 
tetrolet transform, which allows for the acquisition of 
high-frequency and low-frequency coefficients. 
Obtaining the most sparse tetrolet depiction in each block 
is desirable when the tetrolet coefficient is insignificant. 

Step 3: 2*2 blocks are used to reorder the high- 
frequency and low-frequency coefficients. 

Step 4: Accumulation is performed using the high- 


Image 


frequency tetrolet coefficients. 

Step 5: The low-frequency tetrolet coefficients are 
found by going through steps 1-4. 

Output: Decomposed image @ 

The GCV algorithm is used to establish the finest 
denoising threshold value at every scale level in the 
tetrolet field. 

Improved generalized cross validation (IGCV) 

In the wavelet field, the IGCV method is used to 

compute the finest denoising threshold and is also utilized 


in partial differential Here, 


IGCV (6) = arg min MSE(6) if N >0o verifies the 
algorithm. The finest denoising threshold is acquired by 


means of the RGCV method, which is used to evaluate 
MSE. This process is used to maintain the comprehensive 


equations. arg min 


uniqueness of the innovative image. The GCV task is 
used to compute the finest denoising threshold in the 
wavelet field. The GCV principle is used to establish the 
finest denoising threshold to accomplish a_ superior 
denoised image in the tetrolet transform field. 

The input of algorithm is 


y=fte 


(3) 
f — Unknown deterministic structure data. 
€- noise 
y —observation 
W=V+@O (4) 


v=Wf - contains the wavelet co-efficient of original 
data, O=We - Noise co-efficient 

w= Wy 

Here, we select the subsequent GCV task for the 
threshold value: 
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(5) 


N -total number of wavelet co-efficient 

No-—the number of this coefficient that was replaced by 
zero. 

This threshold selection is asymptotically finest for a 
huge quantity of wavelet coefficients. The least Gcv(6) 
also reduces the mean square error task (or risk function) 
R(6) where, 


1 
R(6) =—|lw6 -v]2 (6) 
()= 1 w5—y 
E-denotes the expectation operator in equation (5) 
EGCV (6) = ER(d) +o" (7) 


o° -Denotes the standard deviation 

The higher threshold value contains additional noise. 
The threshold choice process such as the “universal 
threshold” of Donoho and Johnstone is specified as, 


0 = /2log(N)o (8) 


Suppose that X is an innovative image and Y is noisy 
image, then tetrolet decomposes Y whose coefficient 
encompasses low frequency as Y Low, high frequency 
coefficients Y High and allocation of preferred covering 
YC. Afterward, Y Low and YC are retrained and Y High 
is customized by means of threshold which is derived 
from GCV principle. Suppose that if the image is 
decomposed into three stages by means of a tetrolet 
transform, then we identify a matrix m * ny! = y'nign [i,j] 
through the high-frequency coefficients of the row i 
column j in stage |. In the GCV principle, the soft 
threshold is situated as 5 in advance and acquires sj = y |; 
which is the consequence of the threshold. Soft threshold 
task is specified as, 


0 if Vy <o 
S, = ; (9) 
Y, —6 otherwise 


In the tetrolet field, the GCV task is specified as 


2 


1 1 
Fra eed 


GCV(6) = aes 
No 
N 

i is the co-efficient after threshold processing and 


N=m*n, and No=N-S 
S-Number of coefficients that are set to after threshold 


(10) 


processing, 
In this investigational image, the standard deviation of 
Gaussian white noise is 20 and mean is 0. The three-level 
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tetrolet transform is used to decompose the image. The 
GCV (0) of the tetrolet coefficient is used to evaluate 
MSE (6) in dissimilar thresholds. The GCV task is an 
approximate assessment for MSE (6). The threshold is 
equivalent to the least value of GCV (6), which is the best 
denoising threshold. The finest threshold is acquired by 
means of the subsequent equation: 


o,,, =argminGCV(6) 


opt 
Partial differential equation (PDE) 
The Fast Fourier Transform (FFT) process is a 
systematic procedure for arithmetical investigation. It 


(11) 


supplies the estimated investigation for partial differential 

equations (PDE), which is usually premeditated by means 

of the development. It is also known as the basis process, 

which is employed to compute the unidentified 

coefficients. Afterward, we implement the FFT process to 

partial differential equations, which is specified as, 
2 2 

a [St + x) + E(x, y) 

Ot Ox” Oy 


It is used to diminish the number of spatial variables 


(12) 


in anticipation of a two-position boundary or preliminary 
value difficulty, which is resolved by a standard process. 
Here, FFT is used to divide the variables to expand the 
partial differential equation. Due to higher degree, the 
PDE of Fourth order is an aid for diminution and 
division. It is extra flexible in the circumstance and 
division of variables; therefore, the recognition is simple 
and noticeable and diminishes the noise coefficients. The 
identical division is used for the edge circumstances. But, 
it is only relevant for the resolute or limited values 
through the mixture of FFT. Consider © = ©(7,t) the 


field variable (e.g., temperature or concentration), and L 
is a differential operator encompassing additional spatial 
derivatives and a time derivative. Numerous Partial 
differential equations in the beneath 


equation, 


are depicted 


LO= s(r) (13) 

Here S(r) is a task that indicates the dimension 
location. The differential operator is said to be linear i. 

L(a,O, + a,0,) =a,LO, + a,LO, (14) 

Here a; and az are constants ©, (r,t) and ©, (7, t) that 

are tasks which are not essentially an explanation 

L@© = s(r) . We can concern the edge values which must 

satisfy the limit circumstances which are the divided 

values in identical circumstances to investigation. It 


provides the identical occurrence. 
Consider the following PDE, 
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2 2 
: 2 + ¢ . =0 (15) 
Ox Oy 
With the boundary conditions, 
©(0, y) = 0,0d, y) =0 
(0, y) (1, y) 16 


O(x,0) = f(x), O(%,D = f(x) 


The solution can be written in a series expansion as, 


0,9) =¥C,()®, (x) 


n=l 


(17) 


In this equation, 


®,(x) 


®, (x) = V2, sin(na), nl, 2,.... 

Here, we have integrated the Haar wavelet with 
mathematics. The Haar wavelet is a progression of 
task _ that 
constructing a wavelet or source. The Haar wavelet's 


are given by: 


rescaled "square-shaped" consists of 


mother wavelet task of O(t) is specified as, 


1 O<t<l/2 
O(t)=4 -1 1/2<t<1 (18) 
0 otherwise 
Its scaling function g(t) can be described as 
7 1 O<t<l 
a= ( other wise ne 


This least compression proportion accomplishes the 
Harr wavelet coefficients with DWT progression. 
Therefore, the enhanced CR optimizes harr coefficients 
with optimization procedure. In Haar wavelet, the GCV 
algorithm is used to establish the finest denoising 
threshold value at every scale stage. It is an adaptive Haar 
wavelet transformation for image sparse depiction that 
the the 
arrangement and superior energy concentrating capability 
(Mittal, 2023). The foremost problem in the Haar wavelet 
is blocking the object, and GCV will be employed to 
choose the finest denoising threshold in the anticipated 


encompasses uniqueness of geometric 


algorithm. The threshold denoising algorithm is trendy in 
image proposals because it is easy and competent in 
several circumstances. On the other hand, the existing 
denoising threshold choice process is relying on noise 


Statistic possessions like variance. Regrettably, the 
precise arithmetical possessions of noise are normally 
unidentified in definite problems. Therefore, the 


denoising threshold is complicated to evaluate without 
deliberate arithmetical possessions of noise. GCV is 
extensively implemented in image processing, 
establishing an effectual arithmetical process to evaluate 
finest denoising threshold value. The denoised image is 


segmented by means of the Histon-related fuzzy C-means 
algorithm (H-FCM). 
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Histon based FCM 

Image analysis is removing information from an 
image. Image segmentation is a fundamental step in 
automatic pictorial model identification and picture 
investigation, which is the major complicated function in 
image processing. It is the progression of dividing a 
digital image into numerous sections or clusters and is 
employed for establishing objects of concentration and 
limitations such as lines and curves in an image. A group 
of pixels invents every section. The compound image is 
used to discover the function of supplementary image 
processing functions like segmentation, attribute removal, 
machine object identification, etc. Some of the necessities 
of high-quality color image segmentation are as follows: 
A solitary section in a segmented image does not include 
dissimilar colors, and a related section includes an 
identical color that does not encompass additional 
markers. The entire noteworthy pixels contain the 
identical labeled section. The force of the section is 
logically consistent. Color image segmentation is 
constructive in numerous functions. It is feasible to 
recognize the section of concentration and objects in the 
scene from the segmentation consequences, which assist 
the Fuzzy-related 
procedures have an increasing reputation in the area of 


successive image investigation. 
color image segmentation. The anticipated process does 
not necessitate prior information regarding the number of 
objects in an image and computes the threshold value 
mechanically necessary for segmentation. Segmentation 
of compound color image is measured, which unites the 
feature of FCM and A-IFS procedure for to accomplish 
enhanced presentation. This process is used to compute 
association values by means of FCM and A-IFS Histon to 
manipulate the number of clusters. Finally, it utilizes an 
irregularity catalogue to obtain the best threshold values 
for segmentation. The FCM algorithm is implemented to 
discover the association values as depicted in the 
preceding segment. The non-association values and the 
irregularity catalog are established by means of A-IFS 
Histon process. Irregularity catalogue is employed for 
isolating the section. In conclusion, section amalgamation 
is made to accomplish enhanced section related 
segmentation which is reliant on color. 

Generally, digital images are distinguished by means 
of inherent imprecision or ambiguity of grey values. This 
type of uncertainty is maintained by section limitations, 
where it is greatly complicated to encompass an obvious 
separation of boundaries. Clustering procedures are 
derived from FCM, which broadly contracts with 
uncertainty. Intuitionist fuzzy set theory is the trendy 
procedure that is utilized in 


image processing. 
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Additionally, instinctive methods are used to depict the 
inherent uncertainty in images. Intuitionist fuzzy sets are 
anticipated using Atanassov, which is appropriate for 
representing the uncertainty from inaccurate information. 
These are distinct by means of two feature tasks, the 
association task and the non-association task. Association 
task provides a quantity of belongingness of a component 
to the creation of IFS and non-association task provides 
the quantity of non-belongingness. The HFCM depiction 
of images is specified — by means of 


T= 48,4, (8;),vr(ey)} f= be Mf = LeN 


1 =(2,.1(8)).¥;(8))} f= bon Ms Feb 
(20) 
8,- gray values in an image of size MxN, /, (g,;) - 


degree of membership, v, ( 8) degree of non- 


membership. 

The membership value is specified by means of the 
standardized intensity values at every pixel position. The 
anticipated method is diverging from the computation of 
association values. In the anticipated procedure, the FCM 
algorithm is implemented to discover the association 
values rather than the standardized intensity values 
because the association values with the FCM algorithm 
are specified below: 

An image ‘I’ of size pxq with ‘n’ number of clusters 
is considered for analysis. 

(cl,c2c3 ym) L (cl,c2c3) a4) 
I= M O M 
(cl,c2c3), 4) L (cl,c2c3), 59) 


(21) 


The objective function, Jrcy of the image ‘/’ is first 
initialized to zero, 
J rcomo)=0 (22) 
A matrix ‘U’ is used to store the membership 
functions ‘ux,i’ corresponding to each pixel value. 
MaynMen L Meas 
I= M O M 


(Hao Leo L (4) n.0) 


(23) 


The centre ‘0; of each cluster containing the 
membership values of each pixel is given by, 


N 
ee abt od 
i DON on (24) 
eo Axi 
Here v(i=/,2,..., c) is the centre of cluster, 47) 


(x=1,2,..., n, i=1,2,..., c) is association value and Z, = 
(C1, C2, €3) x, (x=1,2,..., n) is the pixel value and ‘m’ is 


fuzzification limitation which is generally obtained as 2 
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for soft clustering. The detachment equation is used to 
discover the distance ‘d’ from the centre of cluster at each 
pixel which is specified by, 

d” (Z,5V;) = (25) 


Using the membership values and Euclidean distance, 


ii 


ty —V,; 


the objective function is calculated by, 
N c in 
pas tied (2,59) =F ney (Y) 


Here, ‘r’ is the number of iterations. The dissimilarity 


(26) 


among the current objective task and _ the preceding 
objective task is specified as, 
Difference = J incu Y)—J ecu ¥ -D (27) 
Suppose that if this dissimilarity is less than the 
threshold value, then the clustering is blocked. This 
threshold value is preferred after creating a widespread 
analysis of a quantity of images to conclude the iteration. 
A threshold of superior value is answered in unacceptable 
clustering for several images. If the dissimilarity exceeds 
the threshold value, the subsequent iteration is established 
from step 4. The innovative association task for next 


iteration is specified as, 
1 


1 (m1) 
d*(z, Ji) 


Xi 1 
x 1 |e 
Al d?(z,_ Ji) 
The average intensity of the pixels is given by, 


Mi (8) = > > my (Gt OG +DAG +k, J+) 29) 


k=l [=-1 
Where, ‘h’ is the filter mask and uw; = UW, ; as found in 


(28) 


equation (29). The hesitancy degree, which lies between 
O and J are given by, 


(8. p))— (8) 
max’; (max, (a (8;,)- i (8:,,))) 


(30) 
The non-membership degree at each pixel location is 
given by, 


m,(g; )=U-“4(8;,;) 


vig )=l-w(s'i-m(s'i)) BI) 

The A-IFS Histon can be defined as, 
¥,(8y) =1- 44 (gif) —2,(8'i/) (32) 
‘g’ lies between 0 to L-J and i = {R, G, B}. ‘6’ 

represents the impulse function. 
d,(m,n) \" 
ionii= os og en) (33) 
o 
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o - Standard deviation of the distance matrix. The total 
distance of all the pixels in a PxQ neighborhood of size 
3x3 is given by, 

d,(m,n) = >) >°dU(m,n)I(p,q)) 


pcPqcQ 


(34) 


The Euclidean distance between two pixels, I (m, n) 
and I (p, q) is given by, 
: (u' (L(m,n))- w' (I(p.q))) + 
d,((m,n).I(p,q)) =| = > (% (l(m,n)) — 3'(I(p,q)))” + 
"(lp g)))” 
(35) 
After calculating the Euclidean distance, a measure 
called roughness index is calculated, which is given by, 


p=1-(f,(g)|+|F(g)) (36) 
Here * F.(g) ” is the histogram of the image. In the 


object section, irregularity catalog is large and it is close 
to concord. In edge section, this value is diminutive and it 
is close to zero. Therefore, an irregularity catalogue is 
valuable in isolating the object section from limitations. 
Finally, the section amalgamation algorithm is used to 
modify the segmented output. After considering a 
quantity of images for investigation, 


Results and Discussions 

To evaluate the performance of the proposed 
technique, the classification is to be performed with the 
aid of the Histon-related fuzzy c means algorithm (H- 
FCM). Here, the performance is to be measured using 
Sensitivity, Specificity, and Accuracy. 


Img4 Imgs 


Img6 


Figure 2. Input Images such as Img], Img2, Img3, 
Img4, Img5 and Img6. 


In the above section, we have included six input 
medical images. 
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Sensitivity 
The measure of sensitivity is the proportion of actual 
positives that are accurately recognized. It relates to the 
ability of tests to recognize positive results. 
Ud os 
(TP + FN) 
where TP stands for true positive and FN stands for 
False Negative 
Specificity 
The measure of specificity is the extent of negatives 


Sensitivity = (37) 


that are properly recognized. It relates to the capacity of 
the test to recognize negative results. 
TN 
(TN + FP) 
Where TN stands for True Negative and FP stands for 
False Positive 


Specificity = (38) 


Accuracy 

The accuracy of the proposed method is the ratio of 
the total number of TP and TN to the total amount of 
data. 


Accuracy = ina ts (39) 
(TN+TP+EFN + FP) 


Discussions 
The simulation results of the proposed and existing 
methods are given below. 


(a) (b) 


(d) (e) (f) 
Figure 3. Img1 (A)Noisy Image, (B) Denoised image, 
(C) Existing K-Means segmented Tumor output, (D) 
Existing K-Means segmented non-Tumor output, (E) 
Proposed Histon based FCM segmented tumor 
output, (F) Proposed Histon based FCM segmented 
non-tumor output. 
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Table 1. Truth table for Experimental outcome. 


: Condition as determined by the 
Experimental 


Standard of Truth 
Outcome 
Positive Negative 
Positive TP FP 
Negative EN TN 


Performance analysis 

In order to evaluate the performance of the proposed 
technique, the classification is to be performed with the 
aid of Histon based Fuzzy C-Means Clustering technique 
here the performance isto be measured by using 
Sensitivity, Specificity and Accuracy. Table 1 shows 
performance measures for our proposed work. 


(a) (b) (c) 


(e) (f) 
Figure 4. Img2 (A)Noisy Image, (B) Denoised image, 
(C) Existing K-Means segmented Tumor output, (D) 
Existing K-Means segmented non-Tumor output, (E) 
Proposed Histon based FCM segmented tumor 
output, (F) Proposed Histon based FCM segmented 
non-tumor output. 


(d) 


Table 2. Sensitivity, Specificity and Accuracy 
Measures for our proposed work. 
Input 


Sensitivit 
Images 7 


Specificity Accuracy 


1 


0.834546337 


0.99527373 


0.990219116 


0.871603623 


0.997486608 


0.993881226 


0.588755856 


0.993146271 


0.981292725 


0.756950953 


1 


0.988128662 


0.924927416 


0.997813861 


0.995 132446 


DW] ny BY] WwW] bd 


0.865736704 


0.985906018 


0.977493286 
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Table 3. PSNR, SSIM and MSE, metrics measures for 
our proposed work. 
Input 


: PSNR MSE SSIM 

images 
1 71.37515734 | 0.004738 | 0.442602877 
2 69.14210699 | 0.007923 | 0.780531334 
3 72.50124437 | 0.003656 | 0.412546045 
4 70.37229725 | 0.005968 | 0.463971721 
5 70.43645289 | 0.005881 | 0.530001656 
6 77.42513233 | 0.001176 | 0.784817566 


(a) (b) 


(e) (f) 

Figure 5. Img3 (A)Noisy Image, (B) Denoised image, 
(C) Existing K-Means segmented Tumor output, (D) 
Existing K-Means segmented non-Tumor output, (E) 


(d) 


Proposed Histon based FCM segmented tumor 
output, (F) Proposed Histon based FCM segmented 
non-tumor output. 


(b) 


(d) (e) (f) 

Figure 6. Img4 (A)Noisy Image, (B) Denoised image, 
(C) Existing K-Means segmented Tumor output, (D) 
Existing K-Means segmented non-Tumor output, (E) 
Proposed Histon based FCM segmented tumor 
output, (F) Proposed Histon based FCM segmented 
non-tumor output. 


Evaluation Measures 


GE Sensitivity 
(ER Specificity 
(RR) Accuracy 


Input Images 
Figure 7. Graphical representation of the proposed 
work in terms of sensitivity, specificity and accuracy. 


From Table 2 and Figure 7, we can understand that 
brain MRI image denoising and segmentation of images 
have good metric values for the BMRI images. Thus, it is 
understood that the medical images can be segmented 
properly by using the proposed work. The proposed 
method achieves the sensitivity, specificity, and accuracy 
measures are 80%, 99%, and 98%. 

80 Mi psnr} 


Evaluation Measure 


1 2 3 4 5 6 
Input images 
Figure 8. Graphical representation of the proposed 
work in terms of PSNR. 


T 


0.4 


Evaluation Measure 


o 9 
NN Ww 


S 


0.0 
1 2 3 4 5 6 


Input images 
Figure 9. Graphical representation of the proposed 
work in terms of SSIM. 
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Evaluation Measure 


2 3 4 5 6 
Input images 
Figure 10. Graphical representation of the proposed 
work in terms of MSE. 
From Table 3 and Figure 8, we can understand that 


brain MRI image denoising and segmentation have good 
metric values for BMRI images. Therefore, medical 
images can be segmented properly using the proposed 
work. Here, we have discussed the PSNR, SSIM, and 
MSE rates and graphically represent the proposed 
measures. 


Table 4. Comparison of proposed and existing 
methods using sensitivity. 


Input Sensitivity for Sensitivity for 
Images _ proposed Histon FCM _ existing k-means 
1 0.834546337 0.867054828 
2 0.871603623 0.937133724 
3 0.588755856 0.585632483 
4 0.756950953 0.664167448 
5 0.924927416 0.915387806 
6 0.865736704 0.87467306 

Comparative analysis 
In this Histon-based image denoising and 


segmentation work is compared with the existing K- 
means clustering technique to prove that the proposed 
work is better for image denoising and segmentation. In 
this section, we have used K-means clustering as the 
existing one. Also, we show in the following Table 6, 
respectively. 

Table 4 and Figure 11 shows that training and testing 
results can be taken for sensitivity. In our proposed 
methodology we have taken input images such as 1, 2, 3, 
45 and 6. The input image 1 attains a sensitivity value of 
0.834546337, and the input image 2 attains a sensitivity 
value of 0.871603623. The input image 3 attained a 
sensitivity value is 0.588755856. The input image 4 
attained a sensitivity value is 0.756950953 and the input 
image 5 attained a sensitivity value of 0.924927416, and 
the final image 6 attained a sensitivity value of 
0.865736704. When we compare the proposed Histon- 
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based FCM results to the existing K-means clustering 
research, our Histon FCM technique will give the best 
result. 


Sensitivity Measures 
o 
Na 


(EG Sensitivity for proposed Histon FCM 
[EE Sensitivity for existing K-means 


1 2 3 4 5 6 
Input Images 
Figure 11. Comparative analysis between proposed 
and existing methods in terms of sensitivity. 


Table 5 and Figure 12 show that training and testing 
results can be taken for specificity. In our proposed 
methodology, we have taken input images such as 1, 2, 3 
4 5 and 6. The input image | attained a specificity value 
of 0.995274 and the input image 2 attained a specificity 
value of 0.997487. The input image 3 attained a 
specificity value is 0.993146. The input image 4 attained 
a specificity value is 1 and the input image 5 attained a 
specificity value of 0.997814, and the final image 6 
attained a specificity value of 0.985906. When we 
compare the proposed Histon-based FCM results to the 
existing K-means clustering, our Histon FCM technique 
will give the best result. 


Table 5. Comparison of Proposed and Existing 
methods using specificity. 


1.00 


0.98 


Specificity Measures 
—) 
ee) 
a 


0.94 


HEE Specificity for Proposed Histon FCM: 
GE Specificity for existing K-means 


0.92 
1 2 3 4 5 6 


Input Images 


Figure 12. Comparative analysis between proposed 
and existing methods in terms of specificity. 


Table 6 and Figure 13 show that training and testing 
results can be taken for accuracy. In our proposed, 
methodology we have taken input images such as 1, 2, 3, 
4 and 5. The input value attained for the accuracy value is 
0.95831. The input image 2 attained an accuracy value is 
0.97168. The input image 3 attained an accuracy value of 
0.969833. The input image 4 attained an accuracy value 
is 0.97168 and the input image 5 attained an accuracy 
value of 0.953827, and the image 6 attained an accuracy 
value is 0.9776. When we compare the proposed Histon- 
based FCM results to the existing K-Means clustering 
method, our Histon FCM technique will give the best 
result. 


Table 6. Comparison of Proposed and Existing 
methods using accuracy. 
Accuracy for 


Accuracy for 


Specificity for Specificity for Input Images proposed Proposed K- 
Input Images _ Proposed Histon existing K- Histon-FCM means 
FCM means 1 0.990219 0.958313 
1 0.995274 0.961276 2 0.993881 0.97023 
2 0.997487 0.971206 3 0.981293 0.969833 
3 0.993146 0.981435 4 0.988129 0.97168 
4 1 0.987471 5 0.995132 0.953827 
5 0.997814 0.955295 6 0.977493 0.9776 
6 0.985906 0.985348 
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0.98 


Aaccuracy Measures 
° 
\o 
ion) 


0.94 


3 Accuracy for proposed Histon-FCM 
(9) Accuracy for Proposed K-means 


1 2 3 4 5 6 


0.92 


Input Images 


Figure 13. Comparative analysis between proposed 
and existing methods in terms of accuracy. 


Conclusion 

In this work, we proposed an effective Histon-related 
image denoising and segmentation algorithm that 
integrates Partial Differential Equation (PDE) and 
Generalized Cross-Validation (GCV) within the Adaptive 
Haar Wavelet Transform framework. The PDE is applied 
to address the issue of blocking artifacts and to estimate 
unknown coefficients, reducing noise more effectively. 
Additionally, GCV is used to determine the optimal 
denoising threshold in the tetrolet transform domain. The 
Fuzzy C-Means (FCM) algorithm was employed to 
discover association values, rather than relying solely on 
standardized intensity values. The experimental results 
show that the proposed method effectively preserves the 
boundaries of digital images. Performance was evaluated 
using Sensitivity, Specificity, Accuracy, PSNR, and 
SSIM metrics. Sensitivity, specificity, and accuracy 
measures for the proposed method were 99%, 98%, and 
98%, respectively, demonstrating significant 
improvement over conventional methods. A comparison 
of the proposed Histon-related FCM approach with the 
existing K-means clustering method further highlighted 
its superiority. Graphical representations of various 
metrics confirm that our proposed method outperforms 
K-means clustering. Therefore, the proposed method 
proves to be a highly effective solution for both image 
denoising and segmentation tasks. In the future, we plan 
to extend this work by exploring more advanced wavelet- 
based techniques and optimization algorithms to further 
enhance the performance of image denoising and 
segmentation. 
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